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The  usual  analysis  of  the  steady  state  angular  motion  of  a  dynamically  stable  spinning 
missile  assumes  a  o^asi-steady  state  calculation  of  a  gravity- induced  trim  angle.  A 
condition  for  the  validity  of  this  quasi-steady  state  assumption  is  derived.  When  this 
condition  is  not  satisfied,  the  gravity-induced  angular  motion  must  be  described  dif¬ 
ferently  for  three  distinct  portions  of  the  trajectory:  the  upleg,  near  apogee,  and 
the  downleg.  The  accuracy  of  this  description  is  checked  by  comparison  with  numerical 
integrations.  Finally  the  influence  of  cubic  static  and  Magnus  moments  on  the  motion 
is  determined  and  a  revised  point  mass  trajectory  model  is  constructed. 
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GRAVITY- INDUCED  ANGULAR  MOTION  OF  A  SPINNING  MISSILE 


ABSTRACT 

The  usual  analysis  of  the  steady  state  angular  motion  of  a 
dynamically  stable  spinning  missile  assumes  a  quasi-steady  state  calcu¬ 
lation  of  a  gravity -induced  trim  angle.  A  condition  for  the  validity 
of  this  quasi-steady  state  assumption  is  derived.  When  this  condition 
is  not  satisfied,  the  gravity- induced  angular  motion  must  be  described 
differently  for  three  distinct  portions  of  the  trajectory:  the  upleg, 
near  apogee,  and  the  dovmleg.  The  accuracy  of  this  description  is 
checked  by  comparison  with  numerical  integrations.  Finally  the  in¬ 
fluence  of  cubic  static  and  Magnus  moments  on  the  motion  is  determined 
and  a  revised  point  r.ass  trajectory  model  is  constructed. 
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1.  INTRODUCTION 


The  linear  angular  motion  of  missiles  can  usually  be  written  as  a 
sum  of  responses  to  various  forcing  functions  and  a  solution  involving 
the  initial  conditions.  For  a  dynamically  stable  missile  the  effect 
of  initial  conditions  quickly  decays  and  the  angular  motion  is  con¬ 
trolled  by  the  forcing  functions,  i.e.  moments  which  do  not  depend  on 
the  missile's  angle  of  attack  or  sideslip  or  their  derivatives.  For  a 
slowly  spinning  missile  the  most  important  such  forcing  function  is  a 
constant  pitch  or  yaw  moment  fixed  on  the  missile  and  caused  by  either 
an  intentional  control  surface  deflection  or  an  unintentional  configu¬ 
rational  asymmetry.  The  response  to  such  a  moment  can  take  on  large 
values  when  the  pitch  rate  is  near  the  roll  rate  and  as  a  result  it 
has  been  studied  by  a  number  of  authors.1-3 

For  a  symmetric  missile  with  a  high  spin  rate,  the  forcing 
function  has  a  magnitude  which  is  proportional  to  the  product  of  the 
spin- to- velocity  ratio  and  the  trajectory  curvature,  and  has  an  axis 
of  rotation  which  is  perpendicular  to  the  plane  of  the  trajectory. 

For  a  constant  amplitude  moment  and  a  linear  static  moment  the  response 
is  a  constant  angle  of  sideslip.  This  trim  sideslip  angle  causes  the 
nose  of  a  spinning  shell  to  always  point  to  the  right  and,  thereby, 
produces  a  right  deflection  of  the  trajectory  which  is  called  drift.4-5 

Since  both  the  spin-to-velocity  ratio  and  the  trajectory  curvature 
increases  to  a  maximum  at  apogee,  a  maximum  gravity- induced  trim  angle 
is  predicted  at  apogee.  This  prediction  assumes  that  a  quasi-steady 
state  calculation  is  appropriate  and  that  the  aerodynamic  moments  are 
linear.  If  either  of  these  conditions  is  not  satisfied,  a  complete 
six- degree- of ^ -freedom  numerical  integration  is  usually  required.  It  is 
the  purpose  of  this  paper  to  present  a  new  simple  approximation  for 


*Referenoea  are  lie  ted  on  page  37. 
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this  gravity-induced  angular  motion  which  is  valid  for  rapidly 
changing  conditions  near  apogee.  The  effect  of  a  non.  near  moment  is 
incorporated  by  use  of  the  quasilinear  assumption  which  has  been  quite 
successful  for  the  analysis  of  the  transient  motion.6'7  Finally  this 
approximation  is  used  to  obtain  a  revised  version  of  a  modified  point 
mass  trajectory. 


II.  EQUATIONS  OF  MOTION 

We  will  make  use  of  two  Cartesian  axis  systems.  The  first  is  an 
earth- fixed  system  with  the  xg-axis  taken  as  the  intersection  of  the 
horizontal  plane  with  the  plane  of  the  trajectory,  the  z  -axis  aligned 
along  the  gravity  vector  and  the  yg-axis  specified  by  the  right-hand 
rule.  The  second  axis  system  has  the  x-axis  along  the  missile-axis  of 
symmetry,  the  z-axis  in  the  plane  of  the  trajectory  pointing  downward 
and  the  y-axis  determined  by  the  right-hand  rule.  For  this  fixed  plane 
axis  system  we  make  use  of  the  complex  angle  of  attack,  £,  which  is 
defined  by  the  equation 


|  =  v-^— —  «»  sin  6  +  i  cos  8  sin  a  (2.1) 

where  v,  w  are  y  and  z  components  of  the  velocity  vector  and  a,  B  are 
the  angles  of  attack  and  sideslip.  The  magnitude  of  £  is  the  sine  of 
the  angle  between  the  missile's  axis  and  the  velocity  vector  and  its 
orientation  determines  the  orientation  of  the  plane  of  this  angle  with 
respect  to  the  horizontal.  For  a  linear  aerodynamics  and  small  geo¬ 
metric  angles  l  must  satisfy  the  equation7 

5"  +  (H  -  iP)|»  -  (M  iPT)g  =  G  (2.2) 

where  the  coefficients  are  defined  in  the  List  of  Symbols. 
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The  plane  trajectory  of  a  particle  acted  on  by  gravity  and  drag 
can  be  described  by  the  equations 


«.  *  -  7  »V*SCD  (r)  (2'3) 

«,  ■  »g  -  j  PV*SC D  (^-)  (2.4) 


Introducing  0^,  the  inclination  of  the  trajectory  with  respect  to  the 
horizontal,  these  equations  can  be  written  in  the  form 


V’ 

V 

=  -  C*  -  g£V-2  sin  0T 

(2.5) 

0rp 

=  -  g&V~2  cos  ©T 

(2.6) 

where 


r*  PS& 
D  *  2m 


C 


D 


Equations  (2.5)  -  (2.6)  can  be  integrated  for  constant  C*  to  give  the 
velocity  as  a  function  of  trajectory  angle. 


V  =  V  sec  eT 
a  T 


t 


tan  0T 


sec  0T  +  In  tan 


(2.7) 


where  V  is  velocity  at  apogee, 
a 
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■MHz] 


~  [3C*  +  4g£V“2  sin  6T]G 


(3.1) 


=  [3C*  +  4g^V'2  sin  8T]PgAV“2  cos  0T 


The  relative  variation  of  G  and  G1  are  indicated  in  Figure  2.  The 
maximum  value  of  G  occurs  after  apogee  due  to  the  action  of  drag.  In 
order  to  obtain  the  angular  response  to  G  we  make  use  of  the  method 
of  variation  of  parameters  for  the  simple  case  of  constant  X^s  and 
(jij's  and  integrate  the  result  by  parts. 


f  =  K  e  1  +  Ke  +  f 
s  2  G 


(3.2) 


=  £ 

G  *g 


r  .  .  (\+iO(s-s)  ,  Cx2+i^)(s-s),  A  a 

Js  [(x2  +  l*i)e  1  1  -  (Xj  ♦  i#J)e  2  JG  (s)ds 

o  (M  +  iPT)  [Xj  -  X2  +  i(*|  -  *2)1 


(3.3) 


For  a  gyroscopically  stabilized  missile  the  fast  rate  is 
usually  much  greater  than  the  slow  rate  •  This  is  especially  true 
near  apogee  where  G*  is  large.  Integrals  of  complex  exponentials  are 
inversely  proportional  to  their  frequencies  and,  thus,  we  can  easily 
neglect  the  first  term  in  the  integral  xn  comparison  with  the  second 
term  and  for  simplicity  we  will  neglect  X^  in  comparison  with  in 
the  multiplying  coefficients  but  not  in  the  exponential  coefficients. 
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(s)ds 


(3.4) 


s 

?  =  „  6  +  i  f 

(*2  +  i^)(s  -  s) 

g  M  J 

6 

0 

L 

where  5  =  . 

g  M 

Although  Equation  (3.4)  is  a  simple  relation  for  the  gravity-induced 
trim  angle  it  is  gravely  limited  by  the  restriction  to  a  constant 
frequency.  For  most  projectiles  the  apogee  value  of  the  gyroscopic 

t 

stability  factor  usually  exceeds  ten  when  G  is  large  enough  to  affect 
Equation  (3.4)  and  a  quite  simple  expression  for  j/  can  be  written  from 
Equation  (2.11). 


* 


1 

2 


M 

P 


(3.5) 


where  s  =  P2/4M. 
g 


Since  P  is  proportional  to  the  spin-to-velocity  ratio  and  the  spin 
normally  decays  quite  slowly  due  to  viscous  damping,  P  can  grow  quite 

f 

rapidly  on  the  upleg  and,  therefore,  the  assumption  of  constant 
is  not  satisfied.  Somewhat  lengthy  algebra  shows  that  a  good  first 
approximation  for  the  effect  of  varying  frequency  on  the  derivation  of 
Equation  (3.4)  is  to  replace  5  by  =  jS  4^  ds  . 

s 

a 


<5  +6  Kre 

g  ga  G 


s)  ♦  i*  (s) 
«  2 


(3.6) 
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where 


KG 


s  -[A  (s  -  sa)  +  i$  (s)] 
|  e 


G’(s)ds 


Equation  (3.6)  clearly  reduces  to  the  quasi-steady  state  relation  when 
O'  can  be  neglected.  Indeed  can  be  neglected  when  the  width  of  the 
humps  of  G  is  large  with  respect  to  the  wavelength  of  $  . 

The  basic  properties  of  Equation  (3.6)  can  be  determined  if  we 
consider  the  very  simple  case  of  zero  drag,  constant  spin  rate  (£  =  0) 
and  no  aerodynamic  damping.  For  zero  drag  Equation  (2.7)  reduces  to 

V  -  Va  sec  6T  (3.7) 


and 


g'  *  4(gtV"2)Ga  cos69T  sin  0T 
Equation  (3.5),  then,  gives  an  approximation  for 


M 

p-  sec  eT  - 
a 


(gt\r2)5-i 

a 


sec  0„ 


(3.8) 


(3.9) 


The  integral  for  now  assumes  a  very  simple  form  for  no  damping 

^2  -i$ 

*G(V  Sg,J  ■  V  /  ■  2  f(«T)d*2  (3-10) 

$2  0 


where  f(8)^  *  4  cos7  6^.  sin  0^  . 
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^2G^ga^  *s  as  a  function  of  <5ga  in  Figure  4.  An  important 

feature  of  this  curve  is  that  K  is  quite  small  for  6  less  than 

2u  ga 

.15  and  thus  we  would  expect  the  quasi-steady  state  results  to  be  good 
when  the  predicted  apogee  steady  state  angle  is  less  than  8°.  When 
the  steady  state  prediction  exceeds  this  value.  Equation  (3.12)  or  the 
more  accurate  Equation  (3.6)  should  be  used.  K„  can  be  identified  as 

di\J 

a  fraction  of  6^a  which  appears  impulsively  at  the  apogee  in  the  slow 

mode  when  the  gravity  forcing  function  is  rapidly  varying  in  a  period 
of  the  slow  mode. 

In  Figure  5  the  combined  pitching  and  yawing  motion  for  a  missile 
with  a  large  maximum  gravity- induced  trim  (2  =  20°)  is  shown.  The 
parameters  used  in  this  exact  integration  are  given  in  Table  I.  From 
Figure  4  we  see  that  is  .76  and,  therefore,  the  amplitude  of  the 

slow  mode  component  of  Equation  (3.12)  is  (.76) (20°)  or  15°. 

This  motion  starts  out  as  an  epicyclic  motion  induced  by  an 
initial  angular  velocity.  During  the  first  three  seconds  the  center 
of  this  epicycle  moves  to  the  right;  then  it  moves  up  as  well  as  con¬ 
tinuing  its  right  motion  until  eight  seconds.  There  then  appears  a 
rough  reversal  of  this  process  until  eleven  seconds  is  reached.  After 
this  point  a  new  epicyclic  motion  is  established  with  a  much  larger 
slow  mode  motion  with  an  amplitude  of  12°.  This  qualitative  behavior 
is  precisely  that  predicted  by  Equation  (3.12)  with  the  near  apogee 
motion  occurring  between  three  and  eleven  seconds.  The  terminal  slow 
node  amplitude  of  12°  is  quite  consistent  with  an  apogee  value  of  15° 
when  the  influence  of  aerodynamic  damping  is  computed. 
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IV.  GRAVITY- INDUCED  TRIM  WITH  DAMPING 

The  effect  of  constant  damping  is  included  in  Equation  (3.6). 
Near  apogee  $2  ®uch  smaller  than  P  and  from  Equation  (2.12)  we  see 

that  a  good  approximation  for  X2  is  -T  which  can  be  constant  for  near 

apogee  flight  For  this  case  Equation  (3.12)  takes  on  the  revised 
form 


5G  =  " 


6  +  6  K_(6  ,  4  )e 

g  ga  G v  ga’ 


■T(s  -  sa) 


(4.1) 


i 

! 


where 


(l|  A  A 

2  T(S  -  s  )  -  i* 

W  |  e  f(ex)d42  -  2»  <  <f>2  <  2n 

-2tt 


On  the  downleg  portion  (  s  >  sQ)  of  the  flight  <j>2  grows  and  H  and 

T  vary  as  the  Mach  number  increases.  During  this  portion  of  flight  G' 
is  quite  small  and  the  integral  in  Equation  i, 3 . 6)  becomes  constant. 

We,  therefore,  assume  the  major  effect  of  G'  is  to  specify  an  initial 
value  of  K2  and  use  Equation  (2.12)  to  predict  the  influence  of 

varying  X2  . 
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r  - 

’G 


i(4  + 

6  +  K  e  2 
g  2 


^.G"1 


(4.2) 


*2  ■  8gaK2G  e*P 


{/ 


(4.3) 


where 


v  ‘*20  ,.-i  ?  T'=  -  sa5  -  i*2(5)  G'  ds 

K2Ge  =ba  |  * 

su 


For  zero  drag  a  simple  expression  for  K2(,  can  be  obtained 


I 


2ir  T(s  -  s  )  -  i|2 


2G  ~  1  ga 


f(0T)  (4-4) 


-2tt 


If  IT/^^ I  <  0.1»  actual  numerical  calculations  show  that  is 


i 

! 

I 


within  .02  of  its  value  for  T  =  0  and,  hence, 

obtain  K„„  as  a  function  of  5 

2G  ga 


Figure  4  can  be  used  to 
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V.  NONLINEAR  ANALYSIS 


The  usual  quasi-linear  analysis6-9  has  been  applied  primarily  to 
the  angular  motions  of  symmetric  missiles  with  no  moment  forcing 
functions.  This  analysis  has  recently  been  extended  to  include  the 
forcing  function  associated  with  slight  configurational  asymmetries.3 
The  latter  treatment  can  be  easily  extended  to  include  gravity-induced 
angular  motion  away  from  apogee.6*10  In  this  section  we  will  outline 
the  appropriate  analysis  and  give  the  results  for  cubic  static  and 
Magnus  moments. 

For  this  case  Equation  (2.2)  becomes 


e"  +  (H  -  IP)  V  -  [M  ♦  M  S2  +  ip  (T  +  T  62)]  i  =  G  (5.1) 

"  ^  o  & 


where 


62  -  m2 

A  solution  of  the  form  of  Equation  (4.2)  is  assumed  and  substituted  in 

i<f>o 

Equation  (5.1).  The  resulting  equation  is  divided  by  K2e  and 

averaged  over  a  distance  which  is  large  with  respect  to  the  wavelength 

of  the  slow  rate  to  yield  quasi-linear  values  of  X  and  . 

2  2 


X 


2 


H$'  -  P  [T  +  T  52  1  + 
r2  0  2  e2J 

2»;  -  p 


(5.2) 


1/2 


-  ^ 


[M  +  M  62  ] 
1  o  2  e2J 


(5.3) 
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where 


<52  =  K2  +  262  +  2K2 
62  2  g  1 

If  the  resulting  equation  is  divided  by  K^e  ,  similar  relations  for 
the  high  frequency  motion  follow.  Finally  the  equation  can  be  averaged 
as  it  is  to  yield  8  quasi -linear  relation  for  the  gravity- induced  trim. 


-  [M 


M  62  + 

2  e3 


iP 


(T  +  T5  . 
o  2  e  3 


)]  5  -  G 


(5.4) 


where 


62  *  $2  +  2K2  +  2K2 

e3  g  2  1 


Since  the  imaginary  part  of  the  coefficient  of  £  is  usually  less  than 

8 

a  quarter  of  the  real  part,  it  affects  the  orientation  of  £  much 

8 

more  than  it  affects  its  magnitude,  6  .  A  simple  equation  for  6 

8  8 

can  be  written. 


G 

6  «  -  (5.5) 

g  M  *  M  (62  +  2K2) 

O  2V  g  2' 

On  the  downleg  Equations  (5.2  -  5.3)  can  be  used  in  Equation  (4.3)  to 
calculate  the  magnitude  of  the  slow  mode  motion  which  has  been  initi- 
ated  by  G  at  the  apogee.  The  orientation  of  the  slow  mode  motion  can 
be  obtained  by  integrating  Equation  (5.3). 

The  nonlinear  analysis  for  near  apogee  motion  is  much  more  diffi¬ 
cult  since  varies  rapidly  during  a  cycle  of  -  An  estimate  of 

the  effect  of  a  cubic  static  moment  can  be  made  for  small  values  of  K„. 
The  steady  state  formulas  for  6^  then  reduce  to  a  cubic  equation. 
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(5.6) 


6  =  - - - 

g  M  +  M  <5* 
o  2  8 


The  slow  frequency,  which  assumes  the  form 


M  +  2M  <$2  1  / 

JL - sec  id-  , 

P  !  T 

a  J 


(5.7) 


varies  in  response  to  the  nonlinearity  as  grows  from  zero  to  6^. 

If  the  nonlinearity  term  in  Equation  (5.7)  is  replaced  by  its  average, 
this  equation  can  be  reduced  to  Equation/  (3.9)  of  the  linear  theory. 


'M  +  M  62  "1 

♦2  *  ~ sec  eT  =  CgfVa2)  6ga  sec  91 
L  a  J 


(5.8) 


Thus,  an  approximation  for  KQ  and  K,^G  when  the  static  moment  is  a 
cubic  function  can  be  made  by  using  Equation  (3.12)  with  6  given  by 
Equation  (5.6)  evaluated  at  the  apogee. 


VI .  A  REVISED  ftJINT  MASS  TRAJECTORY 

For  many  years  ordnance  firing  tables  were  computed  by  use  of  the 
point  mass  Equations  (2.3  -  2.4).  The~e  equations  completely  neglect 
induced  drag  due  to  as  well  as  lateral  drift  caused  by  this  angle. 
The  induced  drag  is  accounted  for  by  adjusting  by  a  form  factor 
which  is  a  function  of  eT  and  is  determined  by  full-range  firing. 

Drift  is  measured  by  full-range  firings  and  numerically  interpolated 
for  firing  table  use. 
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Recently  a  modified  point  mass  analysis  has  been  developed11 
which  includes  the  effects  of  the  steady  state  gravity-induced  trim 
£=-<$.  This  trim  angl^  modified  the  drag  coefficient  and  induced 

o  o 

a  lateral  deflection. 


C  =  C_  +  C_  62 
D  Do  D62  g 


(6.1) 


mye  =  2  p  y2  S  CL  6g 

a  b 


(6.2) 


This  modified  point  mass  trajectory  has  the  advantage  of  retaining  the 
major  trajectory  contribution  of  the  angular  motion  without  requiring 
the  use  of  the  very  small  integration  interval  associated  with  an 
exact  integration  of  Equation  (2.2).  It  is  valid  for  a  dynamically 
stable  missile  and  slowly  varying  G. 

The  theory  of  this  report  can  be  used  to  construct  an  improved 
version  of  the  modified  point  mass  trajectory  which  could  be  called  a 
revised  point  mass  trajectory.  Since  the  motion  near  and  after  apogee 
involves  the  slow  frequency,  an  integration  interval  small  with 
respect  to  the  slow  mode's  period  is  needed.  The  integration  interval 
required  for  Equation  (2.2)  is  small  with  respect  to  the  fast  mode's 
period  and  is,  therefore,  much  smaller  than  that  required  for  the  re¬ 
vised  point  mass  trajectory.  Thu  the  revised  point  rn^ss  trajectory 
requires  much  less  computer  time  than  the  exact  six-degree-of-freedom 
trajectory. 


It  has  been  shown7  that  only  the  average  of  62  need  be  considered 
for  the  drag  force.  If  we  can  neglect  the  effect  of  KG  on  drag  near 

apogee  where  the  total  drag  force  is  small,  only  the  drag  on  the  down- 
leg  need  be  revised. 
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(6.3) 


L 


(62  +  *2) 

g  2J 


The  limitation  imposed  on  the  integration  interval  by  the  slow  mode 
motion  can  now  be  eliminated  if  the  lateral  deflection  due  to  can 

U 

be  neglected.  The  lateral  deflection  due  to  Kg  is  caused  by  an  angle 
which  is  constantly  changing  direction.  The  average  value  of  this  de¬ 
flection  angle  can  be  estimated  by  calculating  the  jump  angle12  for  a 

projectile  performing  coning  motion  of  magnitude  K  and  frequency  <j>'  . 

2  ?» 


This  jump  angle  is  a  right  deflection  angle  of  the  impact  point  with 
respect  to  the  apogee.  The  deflection  angle  with  respect  to  the  gun 
is  one-half  this  angle  and  usually  quite  small.  The  refined  point 
mass  trajectory,  then,  requires  the  same  integration  interval  as  the 
modified  p^int  mass  trajectory.  If  the  effect  of  Kg  on  drag  near 
apogee  is  required,  a  smaller  integration  interval  will  be  required 
for  this  small  portion  of  the  trajectory. 


VII.  COMPARISON  WITH  EXACT  THEORY  EVALUATIONS 

In  order  to  make  a  direct  comparison  with  exact  calculations 

initial  conditions  of  £  =  -6  ,  =  -6'  were  used  with  the  other 

o  g  o  g 
6o  so 

parameters  of  Table  I  to  give  an  angular  motion  without  a  transient  epi¬ 
cycle.  The  total  angle  of  attack  variation  with  time  for  these  condi¬ 
tions  is  shown  in  Figure  6  and  is  compared  with  the  quasi-steady  state 

6  and  the  angular  motion  given  by  Equation  (3.12).  We  see  that  the 
S 

prediction  of  Equation  (3.12)  is  much  better  than  that  of  the  quasi¬ 
steady  state  theory  but  it  does  overestimate  by  35%. 
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The  calculations  based  on  Equation  (3.12)  can  be  considerably 

simplified  if  the  near  apogee  motion  is  approximated  by  a  discontinuous 

jump  at  apogee  from  the  -  <5  motion  before  apogee  to  the  -  6  +  K 

8  8  2 

exp  (i$  )  motion  after  apogee.  This  calculation  which  considers  only 
two  sections  of  the  trajectory  is  also  given  in  Figure  6  and  with  the  ex¬ 
ception  of  a  region  very  close  to  apogee  it  is  seen  to  be  a  good  approx¬ 
imation  to  the  three-section  theory. 

Finally  the  effect  of  damping  is  calculated  through  Equations 
(4.1  -  4.3).  The  two-  and  three-section  calculations  were  repeated 
for  nonzero  damping  and  are  plotted  in  Figure  6.  Here  we  see  that  the 
theory  underestimates  by  about  15%  near  t  =  13  sec.  This  discrep¬ 
ancy,  however,  is  entirely  due  to  calculating  over  a  two-cycle 


interval,  i.e.  one  cycle  on  both  sides  of  apogee. 


was  then  calcu- 

Zu 


lated  over  a  four-cycle  interval  (two  cycles  on  both  sides  of  apogee) 
and  the  result  is  plotted  as  Figure  4.  This  shows  a  difference  of 
about  5%.  The  two-section  calculation  is  repeated  in  Figure  6  using 
the  four-cycle  integration  value  of  K.)(,  and  we  see  the  agreement  with 
the  exact  curve  to  be  quite  good. 
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from  Apogee  for  Several  Values  of  6 


4 CYCLE 


Figure  4.  versus  6  for  X„  =  0.  (For  ! T .2/ 4> ^ l  <  -1  change  in 


» (SEC) 


Figure  6.  Comparison  of  Exact  Values  of  Total  Angle  of  Attack 
with  Various  Theoretical  Approximations 


Table  I.  Parameters  for  Exact 


V 

0 

3 

255  ft/sec 

6 

ga 

= 

eT(0) 

a 

60° 

6 

go 

= 

M 

= 

1.5  x  10-4 

*’ 

?la 

= 

T 

= 

1.5  x  10"4 

4»' 

*2a 

= 

H 

= 

2.9  x  lO’4 

5*  (0) 

s 

P 

o 

= 

.036 

1(0) 

= 

* 

CD 

= 

0 

Integration 

.335  (6  =  20°) 
.071  (Eg  -  1.2°) 
.071 
.002 

-  i  (.5)  rad/sec 
0 
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APPENDIX  A 


DERIVATION  OF  EQUATION  (2.2) 

j 

Most  of  the  relations  given  in  Reference  7  are  for  near  horizontal 
trajectories:  8T  *  0.  These  relations  neglect  the  trajectory  compo¬ 
nent  of  gravity  and  use  an  approximation  for  G.  The  more  exact  equa¬ 
tions  are  given  in  Reference  5  but  in  a  much  different  notation.  It  is 
the.  purpose  of  this  note  to  rederive  the  equations  of  Reference  7  for 
non-horizontal  trajectories.  We  will  refer  frequently  to  equations  in 
Reference  7  and  to  avoid  confusion  will  add  the  letter  R  and  chapter 
number  to  equation  numbers  from  this  reference  to  distinguish  them  from 
equations  in  the  body  of  this  report. 

Equation  (R5.2.2)  is  the  drag  equation  for  a  horizontal  trajectory. 
For  more  general  trajectories,  the  trajectory  component  of  gravity  must 
be  included  and,  thus,  this  equation  must  be  replaced  by  Equation  (2.5) 


-  -  C  J  -  g£\T2  sin  eT  (A-l) 


The  roll  equation  (Equation  (R5.6))  must  now  be  modified  using  the 
more  general  drag  equation: 


♦'*  +  Kp  +  '  -  K6  -  0 


(A-2) 


where  K 
P 


C*  +  cj  +  g*V2  sin  0T 
P 


Preceding  Page  Blank 


39 


k"2  6r  c; 

a  f  £. 


•% 


C 


£ 


6 


f 


+ 


££ 

V 


C„  is  tJ'9  roll  moment  coefficient 
£ 


Foi  a  body  of  revolution,  K  =  0  and 

6 


K  ds, 
P  1 


(A- 3) 


Equations  (R6.4.24-R6.4.25) ,  which  are  the  equations  for  transverse 
motion  with  arbitrary  aerodynamic  force  and  moment,  can  now  be  re¬ 
written  using  the  more  general  drag  equation: 


V 


-  (CD 


♦  gtv  sin  e_  ]  K  -  iy  y 


C-  ♦  i  C~  ♦  ( g-  ♦  i  g~  ) 
Y  Z  \  y  **/ 


£V‘ 


(A- 4) 
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The  non-rolling  coordinate  axes  X  ¥  2  a*e  defined  as  follows:  1)  the 
x-axis  is  aligned  along  the  projectile's  axis  of  symmetry;  2)  the  Y  and 
Z  are  the  orthogonal  axes  of  a  right-handed  Cartesian  set  which  move  so 
that  the  X-component  of  the  angular  velocity  of  the  coordinate  system 
is  zero.  Thus  the  axis  system  pitches  and  yaws  with  the  projectile  but 
has  a  zero  roll  rate. 

A  good  approximation  for  the  aerodynamic  force  is  a  linear  normal 
force. 


Cy  *  1  CZ 


(A-6) 


Equations  (A-4)  and  (A-6)  can  now  be  combined  to  yield  a  revised 
version  of  Equation  (R6.6.1) 


(A- 7) 


The  usual  linear  expression  for  the  aerodynamic  moment  is 


C-  ♦ 

m 


i  C~ 


n 


i 


♦ 


Si 

q 


V  ~ 


a 


(A-  8) 
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Combining  Equations  (A-S)  and  (A- 8)  we  can  obtain 


-  1 


k"2 

Kt 


(A-9) 


Equations  (A- 7)  and  (A-9)  can  now  be  used  to  eliminate  p  and  p'  . 


?,l+(H-~-ip)  l' 


(M  +  i  P  T){  *  G 


CA- 10) 


where 


C* 


c:  -  k 


*  *  <£. )  -  8iv'2 
*  n  1 


sin  0, 


r 


1 

J 


g  =  g;  - 


K2  Si  *  CD ' g£V’2  sin  6t  *  f  * 1  p]s 


l  q 


Gj  =  ^g~  +  i  g;  +  g  Sir.  6T  i'j  £V~2  and 


yCL  =  CN  "  CD 
a  a 


For  small  geometric  angles  (v  =  1,  y'  =  0)  the  definitions  of  H,  M,  and 
T  reduce  to  those  used  in  the  text  of  this  report.  In  order  to  differ¬ 
entiate  ,  we  use  the  identity 


g  =  (kx>  gy.  &l)  +(°>  r)  X  g 


CA— 11) 


Therefore, 


Sy  *  1  8;  ■  i  v 


(A- 12) 


Equations  (A-7)  and  (A-12)  and  some  algebraic  effort  yield  a  good 
linear  approximation  for  G: 


G  =  - 


k‘2  C*  -  C*  .  2g£V“2  sin  eT  +  i  P 

q 


i  P  G 


1 


(A- 13) 
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For  small  amplitude  motion  the  distil*. .ion  between  non-rolling 
coordinates  and  the  fixed  plane  coordinates  vanishes  so  that  tilde 
superscripts  on  £,  ,  y  and  z  can  be  replaced  by  carats  and  G  replaced 

by  G. 


G  = 


i  g  sin 


9 

T  j 


IT 


(A- 14) 


Finally,  the  gravitational  acceleration  perpendicular  to  the  trajectory 
can  be  expressed  as 


g;  =  g  cos 


< 


=  g  cos  8t  cos  a  -  sin  0T  sin 


a 

•m 


g  sin  0 


T 


(A- 15) 


where  =  g  cos  0T 


,  G  =  P 


5  NT 


ig  sin  0T  §  &v 


-2 


(A- 16) 
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APPENDIX  B 


DERIVATION  OF  EQUATION  (3.6) 

Equation  (3.6)  can  be  derived  quite  quickly  by  formally  following 
the  usual  steps  for  the  well-known  method  of  variation  of  parameters 
(see,  for  example,  pages  71-72  of  Reference  7).  The  solution  to  the 
homogeneous  form  of  Equr+ion  (2.2)  with  slowly  varying  coefficients  is7 


B  B. 

-  1  ,  2 
E  =  A  e  +  k  e 
*  1  2 


where  B.  ■  /  X.  ds  +  i  ♦j 


(B-l) 


Aj  are  constants;  and 


Xj  and  *!  satisfy  Equations  (2.11-2.12) 


The  complex  parameters  A^  are  now  made  functions 
particular  integral  for  the  inhomogeneous  term. 
Equation  (B-l),  we  have 


of  s  to  construct  a 
G.  Differentiating 


5  “  A 


B, 

Bj  e  1  +  A 


.  B2 

B  e  2  +  A’ 
2  2  1 


A2  6 


(B-2) 


Preceding  Page  Blank 
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Since  the  differential  equation  is  only  one  condition  on  the  two 
functions,  a  second  condition  can  be  selected.  A  good  choice  is 


B 


+  A, 


CB-3) 


Equation  (B-2)  can  now  be  differentiated  to  yield 


(B-4) 


Equations  (B-l  -  E-4)  are  now  substituted  in  Equation  (2.2)  and  co¬ 
efficients  of  AJ,  A2,  Aj,  and  A;  determined.  Since  Equation  (B-l)  for 
constant  A^  is  a  good  approximate  solution  to  the  homogeneous  part  of 
Equation  (2.2),  the  coefficients  of  Aj  and  A2  must  vanish. 

B  B 

a|  Bj  e  1  +  A’  B’  e  =  G  (B-5) 

Equations  (B-3)  and  (B-5)  can  now  be  solved  for  Aj  and  A ^  and  the 
result  substituted  in  Equation  (B-l)  to  obtain  a  particular  integral 
of  Equation  (2.2).  Indeed,  this  was  the  method  used  to  obtain  Equation 
(3.3)  for  constant  damping  rates  and  frequencies. 

Our  primary  interest  is  in  the  low  frequency  mode  identified  by 
the  2  subscript  and  so  we  solve  Equations  (B-3)  and  (B-5)  for  A^  . 
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Equation  (B-10)  is  equivalent  to  Equation  (3.6) . 
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